THE PARABOLIC-PARABOLIC KELLER-SEGEL SYSTEM WITH 
CRITICAL DIFFUSION AS A GRADIENT FLOW IN R d , d>3 

ADRIEN BLANCHET 1 AND PHILIPPE LAURENCOT 2 

£*\J | Abstract. It is known that, for the parabolic-elliptic Keller-Segcl system with 

critical porous- medium diffusion in dimension M. d , d > 3 (also referred to as the 
quasilinear Smoluchowski-Poisson equation), there is a critical value of the chemo- 
tactic sensitivity (measuring in some sense the strength of the drift term) above 
which there are solutions blowing up in finite time and below which all solutions 
are global in time. This global existence result is shown to remain true for the 
parabolic-parabolic Kcller-Segel system with critical porous-medium type diffusion 

l/"~) ■ in dimension M. d , d > 3, when the chemotactic sensitivity is below the same criti- 

cal value. The solution is constructed by using a minimising scheme involving the 
Kantorovich-Wasserstein metric for the first component and the L 2 -norm for the 

Q, | second component. The cornerstone of the proof is the derivation of additional es- 

^^ ■ timates which relies on a generalisation to a non-monotone functional of a method 

due to Matthes, McCann, & Savare (2009). 
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1. Introduction 
In space dimension 2, the classical parabolic-parabolic Keller-Segel system reads |17| : 



£2: f <v = div[vp-xopVc] 



(t,x) G (0, oo) x 
r d t c = D Ac - a c + fa p , 



in 

en 

en ! where p > and c > are the density of cells and the concentration of chemo- 

attractant, respectively, xo > 0, D , and /3q are positive given constants, and r and 
ao are non-negative given constants. This system is one of the simplest models to 
describe the aggregation of cells by chemotaxis: the diffusion of the amoebae in 
a Petri dish is counterbalanced by the attraction toward higher concentrations of 
chemo-attractant that they themselves emit. This model has been widely studied 
mathematically in the last two decades with a main focus on the so-called parabolic- 
elliptic Keller-Segel system (also referred to as the Patlak-Keller-Segel system or 
the Smoluchowski-Poisson equation in astrophysics) which corresponds to the choice 
r = 0, see [5j HU [151 122] for a review. In particular, a striking feature of the Patlak- 
Keller-Segel system is that, given a non-negative and integrable initial condition po, 
there exists a global solution only if ||po||i < ^XoPo/Do while solutions blow in 
finite time if ||po||i > ^XoPo/Dq. On the one hand, it is worth pointing out that 
this phenomenon has some biological and physical relevance as it corresponds to 
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2 A. BLANCHET AND PH. LAURENQOT 

the formation of aggregates in biology or the gravitational collapse in astrophysics, 
see HU [151 E2] and the references therein. On the other hand, let us also mention 
that such a behaviour is only observed in two space dimensions: indeed, solutions 
do exist globally in one space dimension whatever the value of ||po||i while, in space 
dimension greater or equal to three, there are solutions blowing up in finite time 
emanating from initial data po with arbitrary small mass ||po Hi- 
lt has been shown recently that, in higher space dimensions d > 3, a generalisa- 
tion of this model, which is known in astrophysics as the generalised Smoluchowski- 
Poisson system [12], exhibits a similar threshold phenomenon. It differs from the 
above classical Keller-Segel system by a nonlinear diffusion and reads: 

d t p = div [Vp m - XopVc] , 
(l) { ({,i)e(o,oo)xt d , 

r d t c = DqAc — a c + (3 p , 



where 



m = 2-^G(l,2) 



the parameters Xo > 0) A 5 and (3 are positive given constants, and r and a are 
non-negative given constants. With this choice of m, observe that the porous medium 
diffusion Ap m scales in the same way as the interaction term div [pVc] and that we 
recover m — 1 for d — 2. Several results are now available for the parabolic-elliptic 
version of © corresponding to the choice r = 0, see El ESJ E3 EHJ EHJ EDI EI] and 
it was shown in [HI EI] that there is a critical value M c of the mass ||po||i of the initial 
condition p depending upon d, Xoj A) and (3 such that: 

• if ||po||i < M c , then there is a global solution to (JTJ with initial condition po, 

• given any M > M c , there is at least one initial condition p > with ||po||i = 
M such that the corresponding solution to ([I]) blows up in finite time. 

Though the well-posedness of the parabolic-parabolic Keller-Segel system (CQ) (r > 0) 
has been investigated recently, studies focused on the case where the exponent of the 
diffusion is strictly above 2(d — l)/d [231 ES] and nothing seems to be known when 
m = 2(d — l)/d. The purpose of this paper is not only to show that global existence 
is true in that case for sufficiently small masses ||po||i but also to give a quantitative 
estimate on the smallness condition related to a functional inequality. We actually 
prove global existence as soon as ||po||i < M c , where M c is the already mentioned 
critical mass associated to the parabolic-elliptic Keller-Segel system. 

Before describing more precisely our result and its proof, let us first get rid of 
inessential constants in ([I]) and normalise the initial condition p for p by introducing 
the rescaled functions 

pit, x) = Ru(Tt, Xx) and c(t, x) = Tv(Tt, Xx) , (t, x) E [0, oo) x R f/ , 

with 

n l/(d-2) n d/(d-2) 
D _ n d/(d-2) v _ Up T-A a -nr\ V - R \\n \\ 2 / d 

K — A > A — ~~ u — [Trf - ' 1 — - — ZTU an<1 V ~ Po Po i • 

UPolli UpoII/ 
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Then (u, v) solves 

d t u = div [Vw m — x u ^ v ] 5 
(2) <j (t,x) G (0,oo) x R d , 

rd t v = Av — av + u , 



with 



Xo a ,, ,,2/d . n i a n |i2/d ^ n 

X := -^-PoIIPoIIi >0 and a := d/(d _ 2) IIPolli >0, 
o D Q 



while the initial condition w of u satisfies ||wo||i = 1- Owing to this transformation, 
the smallness condition for global existence on ||po||i is equivalent to a smallness 
condition on x- 

Let us now describe the main result of this paper. We define 

(3) Xc ■= t r-^ , 

(m - 1)Chls 

where Chls is the constant of the variant of the Hardy-Littlewood-Sobolev inequality 
established in [HI Lemma 3.2]: 

h(x) (3^o * h){x) dx 

(4) l7 H ls :=sup<! jRd 2/d : h G (L 1 n L m )(M d ), h ^ \ < oo . 

Ii 

Here 3^o is the Poisson kernel, that is, 

f°° 1 ( \%\ 2 ' 

y (x) := ^ ^^ exp (-_ ) ds = c d \xf ^ 




1 2 — d „ ,- wnd 



where Cd := T(d/2)/(2(d — 2)n d l 2 ). It is shown in [8] that, when r = a = 0, solutions 
to (J2j) exist globally provided that x £ (0, Xc] an d the main purpose of this paper is 
to prove that, when r > and a > 0, solutions to (J2J) also exist globally in time if 
X< Xc- 

Theorem 1 (Global existence). Let r>0,a>0,u bea non-negative function in 
L 1 ^, (1 + |x| 2 ) dx) n L m (l d ) satisfying ||n ||i = 1 and v G H 1 ^). J/x < Xc then 
there exists a weak solution (u,v) to the parabolic-parabolic Keller-Segel system (|2j), 
that is, for all t > and f G C^°(lR d ) ; 

• U6 L oo (0,t;L 1 (R d ,(l + |x| 2 ) dx) n L m (R d )), u m ' 2 G L 2 (0, t; H 1 ^)), 

• u(t)>0, ||«(t)||i = l, 

• ^GL oo (0,t;H 1 (R d ))nL 2 (0,t;H 2 (M' i ))nW 1 ' 2 (0,t;L 2 (M d )) ; v{0) =v , 

and 

[ £ (u(t) -u )dx+ [ [ (V(w m ) - x u Vv) ■ V£ dx ds = , 

jR d Jo JR d 

rd t v — Av + av — u a.e. in (0,t)xR d . 



Note that Theorem [T] is valid without sign restriction on vq. It is however easy 
to check, using the non-negativity of u and the regularity of v, that v o > implies 
v > 0. 
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Remark 2. As for the classical parabolic-parabolic Keller-Segel system in space di- 
mension d = 2, we do not know what happens for x ^ Xc- It is worth mentioning 
that, unlike the parabolic- elliptic case, there might be global solutions to (J2J) for x > Xc 
as recently shown in [I] in the two dimensional case. 

The proof relies on the fact that the parabolic-parabolic Keller-Segel system fl2]) 
can be seen as a gradient flow of the energy 

(5) S a [u,v]:= [ ( l^l" - u{x) v(x) + \ \Vv{x)\ 2 + % v(x) 2 } dx , 

in Vitf^) x L 2 (R d ) endowed with the Kantorovich-Wasserstein metric for the first 
component and the usual L 2 -norm for the second component, where V2(^ d ) is the 
set of probability measures on M. d with finite second moment. Let us mention at this 
point that, since the pioneering works [Td^I^T]. several equations have been interpreted 
as a gradient flow for a Wasserstein distance, see, e.g., [H [2j [TO], [20J [32] , including the 
parabolic-elliptic Keller-Segel system in two space dimensions which actually can be 
seen as a non-local partial differential equation and thus handled as a single equa- 
tion [6j [7], and the literature concerning this issue is steadily growing. In contrast, 
there are only a few examples of systems which can be interpreted as gradient flows 
with respect to some metric involving a Wasserstein distance and, as far as we know, 
a minimising scheme is used to construct solutions to the parabolic-parabolic Keller- 
Segel system (with linear diffusion m = 1) in two space dimensions [TTJ and to a thin 
film approximation of the Muskat problem [IS] . 

To describe more precisely the approach used herein, we introduce the set 

K:= (V 2 nU n ){R d ) xH 1 ^) 

on which the energy £ a is well-defined, see below, and we construct solutions to (J2J) 
by the now classical minimising scheme: given an initial condition (uq,vq) G /C and 
a time step h > 0, we define a sequence {u^ n , %, n ) n >o in K as follows: 

(Uhfi, v h;0 ) = Oo, v ) , 



(6) 
where 



Uh, n +i, Vh,n+i) e Argmin (ttjl , )6X; J'/ l) n[«, v] , n > 



JF h]n [u,v] := — 



+ r \\v-v Kn 



£ a [u,v] 



X 

and W2 is the Kantorovich-Wasserstein distance on T^O^ )■ Several difficulties arise 
in the proof of the well-posedness and convergence of the previous minimising scheme. 
First, as the energy £ a is not displacement convex, standard results from [11132] do not 
apply and even the existence of a minimiser is not clear. Nevertheless, the assump- 
tion x < Xc and a further development of the modified Hardy-Littlewood-Sobolev 
inequality (J3J) allow us to obtain an (L 1 fl L m )(IR d ) x H 1 (IR d ) bound on minimising 
sequences which permits in particular to pass to the limit in the term in £ a [u, v] 
involving the product uv, see Proposition [TJ and prove the existence of a minimiser. 
A similar problem is faced when using this approach to construct global solutions 
to the parabolic-elliptic Keller-Segel system in two space dimensions [61 [7] and is 
solved there with the help of the logarithmic Hardy-Littlewood-Sobolev inequality. 
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A second difficulty stems from the fact that, unlike the parabolic-elliptic Keller-Segel 
system, the Cauchy problem (J2J) cannot be reduced to a single equation. More pre- 
cisely, to obtain the Euler-Lagrange equation satisfied by a minimiser (u, v) of Th,n 
in /C, the parameters h and n being fixed, we consider a "horizontal" perturbation 
for u and a L 2 -perturbation for v defined for 5 G (0, 1) by 

u s = (id + 5(,)#u , vg ■— v + S w , 

where £ G C^°(IR d ;IR d ) and w G CQ°(M. d ) are two smooth test functions, id is the 

identity function of R d , and T#/z is the push-forward measure of the measure \i by 

the map T. Identifying the Euler-Lagrange equation requires to pass to the limit as 

6 ->• in 

Wl(u s ,u h , n )-Wl(u,u htn ) IMg-JIgllg 

_ and _ , 

which can be performed by standard arguments [U [32], but also in 



- / (uv — usvs)(x) dx — J u(x) 



v{x) -v(x + 8((x)) . i .., .. 
i^-J-L - w(x + 6((x)) 



This is where the main difficulty lies since 

wo (id + 5Q — v 

5 ' 



S 



C • Vw in V 



dx 



but u is only in (L 1 fl L m )(IR d ) and m < 2. Therefore, the regularity of (u, v) is not 
enough to pass to the limit in this term and has to be improved. To this end, a 
powerful technique is developed in [20], the main idea being to use as test functions 
in the minimising scheme the solution to a suitably chosen simpler gradient flow 
and analyse the short time behaviour of the outcoming inequality. In the framework 
considered in [20] (see also [7]), the choice of the auxiliary gradient flow naturally 
comes from the existence of another Liapunov functional which is different from the 
energy. Such a nice structure does not seem to be available here but it turns out that 
we are able to somehow extend the scope of this method by finding a simple gradient 
flow (s i—)- (U(s),V(s)) such that s i-> £ a [U(s),V(s)] is the sum of a decreasing 
function and a remainder which can be controlled: for the problem ([2]) under study, 
the choice is the solutions U and V to the initial value problems 

d s U-AU = in(0,oo)xM d , U{0) = u , 

and 

d s V - AV + aV = in (0, oo) x R d , V(0) = v . 

Owing to the properties of (U,V), we have £ a [u, v] < £ a [U(s),V(s)] for all s > 
and this particular choice and the known regularity of (U, V) allow us to pass to the 
limit as s — > in the inequality (S a [U(s), V(s)] — S a [u,v])/s > and finally obtain 
the desired L 2 (M (i )-regularity of u, see Proposition [8] and Corollary [91 This regularity 
allows us to pass to the limit in the Euler-Lagrange equation. This analysis is per- 
formed in Section [2] where we prove the well-posedness of the minimising scheme ([6]) 
and study the properties of the minimisers. Section [3] is dedicated to the proof of the 
convergence of the scheme as h —)■ 0, from which Theorem [1] follows. 
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2. The minimising scheme 

For a > and (u,v) G (L 1 n L m )(M d ) x H 1 ^), we define the functional £ a [u,v] 
by ([5]). First notice that it is well-defined as the integrability properties of u and v 
ensure that uv G L 1 (M d ): indeed, by the Holder and Sobolev inequalities 

(7) ||Hli < \\uhd/(d+2) \\vhd/{d-2) < C ||u||™ /2 \\u\\l /d ||Vv|| 2 < oo . 

We next study the properties of 8 a , a > 0. To this end, let y a be the Bessel kernel 
defined for a > by: 

f°° 1 / Id 2 \ 

the case a = corresponding to the already defined Poisson kernel. For u G L 1 (M d ), 
S a (u) :=y a * u solves 

(8) - AS a (u) + aS a {u) =u in R d 

in the sense of distributions, see [T9| Theorem 6.23]. The Bessel kernel is also referred 
to as the screened Poisson or Yukawa potential in the literature. 

We now show that a modified Hardy-Littlewood-Sobolev inequality is valid for the 
Bessel kernel y a for a > 0: 

Lemma 3 (Hardy-Littlewood-Sobolev inequality for the Bessel kernel). For a > 0, 



h(x) (y a * h)(x) dx 

(9) sup { ^ m : h G (L 1 n L m )(M d ), ft ^ J> = C HLS , 

II A \\m \\h\\ ' 
\\' l \\m H'MIl 

where Chls ^s defined in (jlj). 

Proof. We denote the left-hand side of (jHJ) by /i a . By the definition of y a and 3^o 5 

(10) < y a (x) < y (x) for x6l d . 
Therefore, for h G (L 1 n L m )(M d ), 

ft(x)CV a *ft)(x)dx< / |ft(x)|(;y a *|ft|)(x)dx< / \h(x)\(y *\h\)(x)dx 

JRd JRd 

< (7 \\h\\ m \\h\\ 2 / d 



HLS »|ml"l ' 



whence \i a is finite with \x a < Chls • 

Conversely, given a non- negative function ft G (L 1 fl L' m )(lR d ), we define h\(x) 
\ d h(Xx) for A > 1 and x G M d . Observe that 



'1 1 \ Mr, II _ II I,|| || /, || m \d— 2 MLlim 

tl) flA 1 - ft li ftA \\m — A 



and / ft A (x) (y a * h x )(x) dx = \ d ~ 2 f h(x) (y aX -2 * h)(x) dx . 

JR d JR d 

As A — > oo, y a \- 2 converges pointwisely to ^ for i ^ 0. Moreover, by ffTDl . 
< h(x) (y aX -2 * h)(x) < h(x) (y *h)(x), xeM d . 
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Since h ([Vo * h) G L 1 (R a! ) by the classical Hardy-Littlewood-Sobolev inequality [19 
Theorem 4.3], we infer from Lebesgue's dominated convergence theorem that 

lim / h(x) (y a \-2 * h){x) dx = / h(x) (y * h)(x) dx . 

By (fTTj) . it implies 



h\(x) (y a *h x )(x) dx I h(x)(y aX -2 *h)(x) dx 

/■*<* > lim ^ — - — - — - — -jj- = lim 



\\h>\\™\\h x \\l /d ^°° IHlElHlf 

h(x) (y * h)(x) dx 



H/jllm \\h\\ 2 / d 
\\ u \\m \\ u \\l 

This inequality being valid for all non-negative h G (L 1 fl L m )(K d ), it readily extends 
to all h G (L 1 n L m )(IR d ) so that [i a > Crls , and the proof is complete. □ 

With the above notations, we have an alternative formula for S a for a > 0. 

Lemma 4. Let a > 0, u G (L 1 n L m )(R d ), and v G H 1 (R d ). Then 

£ a [u,v] = £ a [u,S a (u)} + -\\V (v - S a (u)) \\l + -\\v- S a (u)\\l. 

Proof. We proceed as in [9] and first claim that S a (u) G H 1 (M a! ) so that S a [u, S a (u)] 
is well defined. Indeed, since m G (1, d/2) and u G (L 1 fl L m )(M. d ), we infer from [T9| 
Theorem 6.23 (v)] that S a (u) G V(R d ) for all r G [l,dm/(d - 2m)]. In addition, 
using a rescaling with respect to the parameter a, we deduce from [251 Chapter V, 
§ 3.3, Theorem 3] that, for all p G (l,m], there is C = C(d,p) such that 

(12) all^^Hp + v^llV^^llp+IID^^IIp^CCp) \\u\\ p . 

In particular, since 2d/(d+ 2) G (1, m), we have S a (u) G W 2,M/(d+2) (IR d ) and the 
Sobolev embedding guarantees that VS a (u) G L 2 (IR d ). Also, dmjid — 2m) > 2 and 
we thus have S a (u) G L 2 (M a! ), thereby completing the proof of the claim. 
We now compute 

£ a [u,v]-£ a [u,S a (u)]= / { -u[v- S a (u)} + - \V[v - S a (u)}\ 2 

+ V[v - S»] • VS Q (u) + | k - S»| 2 
+ a [u — S^u)] ^(m) f dx 

= 1 -\\V[v-S Q (u)]\\l + ^ 



^ \\V\v - S a (u)]\\l + ~ ||u - S a (u) n , 



+ I [v — S a (u)] [—AS a (u) + aS a (u) — u)] dx 

JR d 

which gives the stated result by using the Bessel equation (jSj). D 

A lower bound on £ a follows from Lemmas [3] and HJ 
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Lemma 5 (Lower bound on £ a ). Let Xc be defined in fl3]). If a > and (u,v) G 
(L 1 n L m )(R d ) x H^R*), then 

/■\ c r i \ ^ HLS ( II n 2 /<A II \\m 

00 £a[u,v]>— [Xc-X\\u\\l Mm ■ 



2X 
Moreover there exists C\ > depending only on d such that 

(ii) II Vv||| + al|v||l < 4 f„[w, v] -h C7i \\u\\T \\u\\Z ■ 

Proof. Owing to the property (J8j) of the Bessel potential and the modified Hardy- 
Littlewood-Sobolev inequality for the Bessel kernel ([H]), we have: 

Hull™ 1 f 

S a [u, S a (u)] = ." "" - - / u S a {u) dx 
X(m - 1) 2 J Rd 

\\u\'~ 

> 



X(m- 1) 
. Chls ( 



<--HLS 
2 


\\u\\ 


m 


U 


.2/rf 

1 


-x N 


,2/rf' 


) 


l«|| 


in 
m ' 



2X 

The estimate (Ej) then readily follows from Lemma H] and the above inequality 
Next, by (JTj) and Young's inequality, we have 

|| Vf || 2 + « || ^|| 2 < 2£ a [u, u] + 2 1| wu ||i 

<2£ Q [u,i;] + C ||n||™/ 2 ||u||; /d ||W|| 2 

<2S a [u 1 v} + ^^ + C\\u\\ 2 1 /d \\u\\ m 



2 



jit 



<2R,K»l+ " V '" i ^'' IHIi +0 11.11^ N|g, 

which gives the stated result (Jn|). D 

Though Lemma H] is not true for a = since Sq(u) $_ H 1 (IR d ), it turns out that 
Lemma [5] is still valid in that case as we show now. 

Lemma 6 (Lower bound on £ ). Let Xc be defined in ([3]). If(u,v) G (L 1 flL m )(]R d ) x 
H 1 (IR d ) ; then the statement of Lemma\^ remains true for a = 0. 

Proof. Since we cannot use (u,S (u)) in £ , we use an approximation argument: 
recalling that, given e G (0, 1), S £ {u) G H 1 (R d ), we may compute 

£ [u,v]-£ [u,S £ (u)} = l\\V(v-S E (u))\\l-e [ S £ (u) (v - S £ (u)) dx . 



2 

We infer from Holder's, Sobolev's and Young's inequalities that 
£ [u,v] -£ [u,S £ (u)] 

> - \\V(v-S e (u))\\l-e \\S £ (u)\\ 2d/(d+2) || (v - S £ (u))\\ 2d/{d _ 2) 

> \ W{v-S e {u))f 2 -Ce \\S £ (u)\\ 2d/{d+2) \\V(v-S £ (u))\\ 2 

> —C e \\S £ {u)\\ 2d /( d+2 j , 



S e (u) (u + e S £ {u)) dx 
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whence 

(13) S [u,v] > S [u,S e (u)] - C e 2 \\S e (u)\\ 2 2d/{d+2) . 

It next follows from the modified Hardy-Littlewood-Sobolev inequality ([9]) that 

\\ll\\ m 1 f 

£o[u, S £ (u)} = l lu - - / S e (u) (2u + AS e (u)) dx 

X{m — i) i j R d 

X(m-l) 2 

> \\ U \\m _ CHLS I, |, m || „2/d _ £ || o / A..2 

- X (m-1) 2 HLIMIi 2 IIMuJlla 

. C'hLS / I, ||2/<A || ,, m £ n c / n|| 2 

-~2 — ^c-xlHli J ML- 2 H&wlU- 

Combining (ITBl and the above inequality gives 

(14) £oM>^- (xc-xll«li; /d ) IW-Ge a ||5 , .(tt)||^ 2) -|||5'.(«)||2. 

We are then left to study the behaviour of the last two terms of the right-hand side 
of ( TT4"j) as e — > 0. To this end, we note that S^w) = 3^o * ( M ~ £ S £ (u)) and it follows 
from the regularity of solutions to the Poisson equation (see, e.g., [T9| Theorem 10.2]) 
and the estimates ffT2l) for the Bessel potential that, for p G (1, m], 

(15) \\S £ (u)\\ dp/{d ^ 2 p) < C(p) \\u-e S £ (u)\\ p < C(p) \\u\\ p . 

We set p m := dm/(d — 2m) > 1 and consider p G (1, 2d/(d + 2)) to be determined 
later on. Since 2 G (j>o,Pm) and 2d/(d + 2) G (po,Pm), we have 

A^-LJ5L e(0 ,i), /.== M M -y + 2)a> , 6(0.1). 

Pm-Po (d + 2)(p m -p ) 

and we infer from Holder's inequality, ( 1121) . and (114J1 that 

e||5' £ ( U )|| a < £ ||S' £ ( W )|t ||^(«)||£- A)P0 

(16) < C7(p ) ^(a-i)^ || u ||£™ || m ||(i-a)po 



I PO 



and 

e 2 



s'.WIImaw) < ^ 2 ll^(^)llii +2Wm/d ll^(«)llg +2)(1 -^ 0/d 

(2d+{d+2)(ix-l)po)/d u u(d+2)tip m /d || ||(d+2)(l-At)po/d 



;>u 



(17) < C(po) £i 2< *+i d +w-i;po)/<* |U||(d+2) AiPm /d n 

Since 

1 , C\ i\ Pm+P0-PmP0 _^ 1 ^ n 

i + (A - i)po = — > > , 

Pm ~ PO PO^ 1 Pm ~ 1 

and 

9 , x//x9 w n p m (2cj-(d + 2)p ) (d-2) Pm ^ 
2d+ (d + 2)(fx - l)p = — > — >0, 

Pm - PO PO^ 1 Pm - 1 

we can find p G (l,2d/(d+ 2)) sufficiently close to 1 such that 1 + (A — l)p > 
and 2d+ (d+ 2)(/i — l)p > 0. Owing to ([TBI and (fT7j) . this property readily implies 
that the last two terms of the right-hand side of (|14p converges to zero as e — > 0. 
Consequently, letting e — > in ( TT4l) gives that the estimate (El) of Lemma [5] is also 
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true for a = 0. The proof of the estimate for ||Vw||2 is then the same as that of 
Lemma [5] (jnj) . □ 

After this preparation, we are in a position to define and study the minimisation 
scheme. Set 

K:= (7> 2 nL m )(R d ) xH 1 ^), 
where T > 2(J^ d ) is the set of probability measures on IR d with finite second moment. 
For a > 0, h G (0, 1) and (u ,v ) G /C, we define the functional 

1 /W|(«,«o) , l|2 



F h [u,v] := — I ^t\\v-v \\ 2 \ +£ a [u,v] , 

where W2 is the Kantorovich-Wasserstein distance on P 2 (K d ). Our aim is to minimise 
JFh on /C. 

2.1. Existence and properties of minimisers. 

Proposition 7 (Existence of minimisers). Given \ G (0, x c ), (u ,t>o) G /C, and 
/i G (0,1), i/iere exists at least a minimiser (u,v) G /C of Fh in /C. Moreover, any 
minimiser (u,v) of Th in K, satisfies 

(18) T h \u, v] < F h [u , v ] = £ a [u , v ] . 

Proof. We proceed in two steps. 

Step 1: Estimates. Consider (u, v) G /C. The lower bound on the free energy, see 

Lemma [5] (i) , ensures that 

(19) F h [u, v]> — ( 2K 0) + r\\v - v \\l )+ rj\\u\ 



m 5 



2h V X 

where 77 := Chls (Xc ~ x)/^X > 0. Therefore, Th is non- negative in K and we may 
define 

(20) u := inf F h [u,v] > . 

(u,v)eK 

Let (uk,Vk)k>i be a minimising sequence with a; < ^[iifc, ffc] < w + 1/A; for fc > 1. 
Recalling that 

m(S ,fi) < m(5 0l u ) + W 2 ( Mo ,/i) , H G P 2 (M d ) , 

where 5q denotes the Dirac mass in M. d centred at x = 0, we deduce from ( Tl9"l) that, 
for fc > 1, 



|x| Uk(x) dx < 2 / |x| Mo(^) dx + 2W 2 («fc, Mo) 
< 2 / |x| 2 u (x) dx + 4hxFh[uk, v k] 

JR d 

(21) <2 \x\ 2 u (x) dx + Ah X (uJ + I) ■ 



It also follows from $gfy and flUD that, for k > 1, 



„ Fh[uk,Vk] , 2/iJ) l [n fc ,f fc ] ^ (1 t 2h\ 



(22) ||tt fc ||™ + || V|k - v U < ^i^l^i + — "t"«."M < - + - (w + 1) 

77 r \rj t J 
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Furthermore, we infer from Lemma [5] (ii) that 
(23) \\Vv k \\i<AE a [u k ,v k }+C 1 |KC<C, k>\. 



11 



Step 2: Passing to the limit. Owing to (ED), (J22D, and (ESj), it follows from the 
Dunford-Pettis theorem that there are (woo,«oo) G (L 1 H L TO )(R d ) x H 1 (E d ) and a 
subsequence of (u k ,Vk)k, which is not relabelled, such that 



(24) 
(25) 



i'k 

Vk 



II r 



in (L 1 n L r 



inH 1 



Since Uk is a probability measure for all k > 1, the convergence (1241) guarantees 
that Uoo is a probability measure. Next, a classical truncation argument and ( 1211) 
imply that the second moment of Mqo is finite. Therefore, Mqo belongs to T^C^ ) and 
(tioc^oo) e /C. 

Next, weak lower semicontinuity arguments and the properties of the Kantorovich- 
Wasserstein distance allow us to deduce from (1241) and (1251) that 



(26) -( +T\\ Veo 



VoM + 



i«ooii; 



|VVoo|ll 



a m 



2 
oo||2 



X(m-l) 2 2 

< liminf I J^K,^] + / u k (x) v k (x) dx 

k^oo \ J Rd 



We are then left with passing to the limit in the last term of the right-hand side 
of (126]) . For that purpose, given n > 1, we fix a truncation function 9 n e C^° r 
satisfying < 6 n (x) < 1 for all lei" 1 and 



n 



9n(x) 



if \x\ < n 
if \x\ > In . 



By (|23|) . {9 n Vk)k>\ is bounded in 14^(5(0,277,)) for each n > 1 and, using a diagonal 
process, we can extract a further subsequence, not relabelled, of (v k)k>i such that 

(27) 6 n v k -> # n Woo in L p (M d ) for any p G [2, 2d/(d - 2)) and n > 1. 

It next follows from the Holder and Sobolev inequalities that 

(d+2)/2d 

■ : 

< IK ; fc||2d/((i-2) 



6> n ) u k v k dx 



(« 



U{M>n} 

<C||Vi; fe || 2 ||n fc C/ 2 
^l|V^|| 2 |K||™ /2 



«f (d+2) dx 



life dx 



L/{H>n} 



i/d 



< 



71 



\x\ 2 Uk dx 



l /d c 

< 



n 



2/d 



the last inequality being a consequence of (|2T|) and (122]) . Similarly, since (moo,w c 
/C, we also have 



(29) 



^n) "oo ^oo U-X 



< 



c 

2/d ' 



'/? 
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Thanks to (122]). (125]). and (1291). we have for all n > 1 
(w fc v k - Moo Woo) dx / 9 n (v k - v^) u k dx 

(1 - 9 n )u k v k dx 

< C \\9 n (v k — foo)||m/(m-l) + 

+ Cn- 2 ' d . 



n (Uk - Moo) Woo dx 



(1 - 6'„)n 00 t; 00 dx 



n (lift - Moo) Woo dx 



Since 2 < m/(m - 1) < 2d/(d - 2) and OnV^ G L m /( m " 1 )(M d ), the first two terms 
of the right-hand side of the above inequality converge to zero as fc — >■ oo by 
and (12TI) . Thus, for all n > 1, 



lim sup 

fc— >oo 



(WfcWfc - Moo ^oo) dx 



< Cn~ 2 / d 



Letting n — > oo gives 



(30) 



lim / Mfc(x) v k (x) dx 

fc— >oo 



u 00 (x)v 00 (x) dx 



Combining the inequalities (|26|) . (|30|) obtained previously, and the minimising prop- 
erty of (u k ,v k ), we end up with J-'hiuoo, Woo] < w. Since (moo,Woo) belongs to /C, we 
actually have J-'/Jwoo, Woo] = w by (1201) and (woo, foe) is thus a minimiser of J^ in /C. 
Finally, given a minimiser (u, u) of J-jj in /C, the inequality (IT8"1) is a straightforward 
consequence of the minimising property of (u, v). □ 

We next improve the regularity of minimisers of T^ in /C by adapting an argument 
developed in [2U|. The key argument in [2U] is the existence of an additional Liapunov 
functional for the problem under study which is associated to another "simpler" 
gradient flow, the solutions to this gradient flow being then used as test functions. 
Though there does not seem to be such a structure hidden in the parabolic-parabolic 
Keller-Segel system ([2]), we nevertheless find a simple evolution system, the solutions 
of which we use as test functions and thereby obtain additional regularity for the 
minimisers. 

Proposition 8 (Further regularity of the minimisers). Let \ £ (0, Xc)> ( u o 5 v o) £ £■, 
and h G (0, 1). Consider (u, v) G K, be a minimiser of Th inJC. Then, u m l 2 G H 1 (IR d ), 
An — an + u G L 2 (M a! ) and there exists C 2 > depending only on d, x, a , an d T su °h 
that 



(31) 



mx 



\W(u m / 2 )\\ 2 2 + \\Av-av + u\\l<2A h (0) + C 2 (£ a [u ,vo}+£ a [uo,v } 1 ^ m -^) 



where 



(32) A h (0) := V^A_2M + 1. (|| VvQ g + a\\v \\ 2 2 - \\Vv\\l - a \\v\\ 2 2 ) 



THE PARABOLIC-PARABOLIC KELLER-SEGEL SYSTEM AS A GRADIENT FLOW 13 

and % is Boltzmann's entropy defined by 

(33) %\w] :— / w(x)\og(w(x)) dx . 



Recall that, if w e {V 2 n L m )(M d ), then wlogw e L 1 (M d ) and there is C 3 > 
depending only on d such that 



w \logw\ dx < C 3 (l + ||w||m) + / w(x) (1 + \x\ 2 ) dx 

JR d 

(34) H[w] > -C 3 - w{x) {1 + \x\ 2 ) dx , 



see, e.g., [THJ p. 329]. Therefore, *H[v,o] and K[w] are well-defined and thus Ah(0) as 
well. 

Proof of Proposition^ Let (u,v) be a minimiser of 7-^ in /C. We introduce the 
solutions U and V to the initial value problems 

(35) d t U - AU = in (0, oo) x M d , 17(0) = « in R d , 
and 

(36) d t V -AV + aV = in (0, oo) x R d , 1/(0) = t) in M d . 
Classical results ensure that (U(t), V(t)) belongs to /C for all t > and therefore 

(37) F h [u,v]<F h [U(t),V(t)], t>0. 

Unlike in [20], the free energy S a is not a Liapunov functional for (U, V). Nevertheless, 
and this is actually the cornerstone of the proof, the time derivative of S a along the 
flow of fl35|) - fl36|) is the sum of a negative term and a remainder term which can be 
controlled. The claimed additional regularity on (u, v) results from the computation 
of the time evolution of .7-^ [77, V] and its proof requires several steps: we first use ( 1331 ). 
( 1331) . and their gradient flow structures to compute the time evolution of J-^U, V]. 
The second step is devoted to control a remainder term arising in this computation. 
Once this is done, weak convergence arguments and the lower bound on S a , see 
Lemma [5] and Lemma El are used in the last step to obtain (1311) . 
Step 1. 
• It follows from (133)) . (1331) . and integration by parts that 

-^-£ a [U,V] = -— [ U m ~ 2 \VU\ 2 dx- [ U (AV-aV + aV) dx 
- f U (AV - aV) dx- f (AV - aV) 2 dx 

JR d JR d 

where 



and 



V(t) := — ||V {U m ' 2 (t)) ||i + ||(AV -aV + U)(t)\\l , * > 
mx 

n{t) := \\U{t)\\ 2 2 -a I (UV)(t,x)dx, t>0. 

JR d 
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After integration we obtain 

(38) £ a [U(t), V{t)} - £ a [u, v}<- V{s) ds + / K(s) ds , t > . 

Jo Jo 



• We next recall that the linear heat equation (135]) can be interpreted as the gradient 
flow of the functional H defined in (13B1 for the Kantorovich-Wasserstein distance 
W2 in P 2 (^ )> see [161 EI]- Therefore, owing to the differentiability properties of 
the Kantorovich-Wasserstein distance (see, e.g., [U Corollary 10.2.7] and [32], Theo- 
rem 8.13]), it follows from [U Theorem 11.1.4] that 

l^2(U(t),u ) < H[u ] - H[U(t)} , t > . 
Integrating with respect to time we obtain 

^ [W 2 2 (U(t),u ) - W 2 (u,u )] < J (H[u }-H[U(s)}) ds . 

By the monotonicity of s i-> T-t[U(s)] we deduce that, for all £ > 0, 

(39) 1 [W 2 (£/(£), wo) - W 2 (u, n )] < t (U[u ] - H[U(t)}) . 

• Furthermore, it readily follows from ( |36l and Young's inequality that 

o-dl^-^o||i = - / [V(V-v )-VV + aV (V-Vo)} dx 
z at j Rd 

<-\\VV\\i + ±(\\VV\\i + \\Vv \\l)-a\\V\\l 

+ f (l|V||2+IMIa) 

<-\(\^V\\l + a\\V\\l)+l(\\Vv \\l + a\\v Q \\l) . 

Integrating with respect to time and using the monotonicity of s *— > ||W(s)||| and 
s i-> ||V(s)||| we end up with 

(40) \\V(t) - v \\l - \\v - v \\l < t (\\Vv \\l + a \\v \\l - \\VV(t)\\l ~ « ||^(*)||1) 

for all t > 0. 

• Combining the above estimates (137|) . (138|) . (1591 and (1401 gives, for £ > 0, 

0<F h [U(t),V(t)]-F h [u,v] 

< -^- (H[uo] - n[U(t)}) - [ V(s) ds + / K{s) ds 
^X Jo Jo 

+ ^ (||V«b||2 + a M^olll - || W^H 2 . - a ||V(£)|| 2 ) , 



which also reads 



1 /"* 1 '* 



(41) 7 / D(s) ds < A h (t) + - / ft(s) ds , t > 

'0 Uo 
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with 

Mt) ■= n[U ° ] ~^ [m] + Yh (II V^olll + a \\v \\l - \\VV(t)\\l - a \\V(t)\\l) • 

Owing to the continuity and regularity properties of U and V, the function Ah has 
a limit as t — > and 

(42) lim A h (t) = A h (0) , 

t-s-o 

the constant Ah(0) being defined in ( 1321) . 

Step 2. We now estimate TZ. Since 2 G (m, 2/(m — 1)), the Holder and Sobolev 
inequalities ensure that, given iy G L m (M a! ) such that |w;| m / 2 G H 1 (R d ), the function 
w belongs to L 2 (IR d ) and 

(43) \\w\\l < \\w\\ m |MU /(m _i) < C \\w\\ m ||V(|«;r /2 )||2 /m . 
Combining the above estimate with Young's inequality gives 

\\u\\\ < — \\v(u m/2 )\\l + c ||c/||^ /(m-1) , 

mx 
and thus 

(44) \\U\\l<^ + C\\U\\Z l{m - l) . 
Moreover, we infer from (J7J) and Young's inequality that 

(45) a\\UV\\ x < aC\\U\\Z 12 \\U\\\ /d || W|| 2 < i^S + o?C\\U\\2 ■ 

We now infer from ( 14*4|) . (|45l) . and the time monotonicity of s 4 || W^(s) ||| and 
s^||C/( S )||£that 

^V( S )d S <i^( S )d S + i||V^|| 2 + C(||n||-+||n||^- 1 )) . 

By OH]), Lemma 0and Lemma M we finally obtain 

]fv(s)ds < 2A h (t) + \\Vv\\l + C(\\u\C+\\u\C« m -V) 

(46) < 2A h (t) + C(E a [u,v]+S a [u,v} 1 / ( - m -V) . 

Step 3. We are now left with passing to the limit as t — > in the left-hand side 
of (146]) . To this end, define first 

1 /"' 
V x (t,x):=- l U m/2 (s,x)ds, (t,x) G [0,1] xR d . 
t Jo 

Since m G (1, 2) and U G C([0, 1]; L m (R d )), the function V x belongs to C([0, 1]; L 2 (M d )) 
and T> 1 (t) converges to u m ^ 2 in L 2 (IR d ) as t — > 0. Combining this property with (142)) 
and (|46p . we realise that (VDi(t))te(o,i) is bounded in L 2 (R d ;R d ) and converges to- 
wards V (w m/2 ) in H- 1 (M d ;M d ). Therefore, V (w m/2 ) belongs to L 2 (M d ;M d ) and 

(47) ||V(n m / 2 )||' < liminf l|VX>i(t)l|2 . 

II ||2 t _^ nz 
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Similarly, defining 

V 2 (t,x):=- ( A / V(s,x) ds-a V(s,x)ds + U(s,x) ds 
t V Jo Jo Jo 

for (t,x) G [0, 1] x R d , we infer from (JI2D and flU) that (P 2 (*))te(o,i) is bounded in 
L 2 (IR d ) while the continuity properties of U and K with respect to time guarantee 
that T> 2 {t) converges towards Av — av + u in H~ 1 (M d ) as t — >■ 0. Consequently, 

Aw - aw + u G L 2 (IR d ) and 

(48) ||Av - av + u||2 < liminf ||X> 2 (t)||2 • 

Thanks to (02)), flU), and (JMD, we may let £ -> in Q3SD and obtain the stated 
result. □ 

Corollary 9. Let x £ (0,Xc), (uo^o) £)C, /i G (0,1), arad consider a minimiser 
(u,v) ofT h in K. Then u G L 2 (M d ). 

Proof. Since u belongs to L m (R d ) by the definition of /C, Corollary follows at once 
from Proposition 181 and (l4"3"j) . □ 

2.2. The Euler-Lagrange equation. 

Lemma 10 (Euler-Lagrange equation). Consider x £ (0,Xc), (wo,fo) G ^, and 
/i G (0, 1). If (u,v) G /C is a minimiser of Th in /C, then 

Wi(u,u ) 



(49) 



/ [£(u -u ) + hV£- (V (w m ) -p Vv)] dx 



< lien 



W 2 -' 



V — Vq 

(50) r Av + av — u = a.e. in 

h 



for any £ in C^°(M. d ) and 

(50) r — f 
In addition, V (u m ) ~xuVv e L 2 ™^ 2 ™" 1 ) (M d ) and satisfies 

(51) & ||V(u™)-x« Vw|| 2m/(2m _ 1} <C 4 W 2 (u,u ) ||V(^ 2 )|| 2 /m 
for some positive constant C^ depending only on d, a, x an d T - 

Proof. Pick two smooth test functions ( G C~(R d ;R d ) and w G Co > (M d ) and define 
Ts := id + 5 ( and 

u s = T s #u , v§ ■= v + 5 w 
for 5 G (0, 1), where id is the identity function of lR d and T#fi denotes the image 
measure or push-forward measure of the measure [i by the map T. Notice that there 
is 8{ small enough such that Ts is a C^-diffeomorphism from M. d onto M. d for all 
5 G (0, o" c ) and Det (VT S ) = Det (1 + 5 VC) > 0. 

• By a standard computation, see [321 Theorem 5.30] for instance, 

||„. Ilm IL. ||in p 

(52) lim " 5lm ' " m = - / Tr(VC(x)) u m (x) dx . 
s^o (m - l)d J Rd 

• It is also standard, see [321 Theorem 8.13] that 

f K o\ v Wi(u S ,U ) -Wi(u,U ) f ( X7 ( \\ f X7 t \ ( \A 

(53) hm — = — / [x — \/<f(x)) ■ C° v<p{x) Uo{x) dx , 

(5-5.0 26 . md 
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where V</3 is the optimal map pushing uq onto u, that is, u = Vf^uo and 
Wliu.uo) = / |x — Vy?(x)| 2 Uq(x) dx 



= inf < / |x — T(x)| u (x) dx : T#u = u 

{jR d 

• Another classical computation ensures that 

(54) lim— [HVW5II2 + a|Kll2 _ \\^ v \\l - a \\ v \\l] = / (Vv-Vw + avw) dx 

• Finally, by the definition of the push- forward measure, 
u v — us vg) (x) dx = / u(x)[v(x) — v(Tg(x)) — 5w(Ts(x))} dx . 

JR d 

Since u is bounded in L 2 (M (i ) by Corollary [9] and 

C-Vv in L 2 (M d ), woTg ->• to in L 2 (M d ), 



«°?>-V vA.TT.. ;„r2^ 



S 
we conclude that 



(55) lim - / [u v — ug vg] (x) dx = — / (u( -Wv + uw) dx . 

• Since (ug,Vg) belongs to /C, we infer from the above estimates (152]) . (1531) . (E 
and fl55l) that 



< lim - (J>[u,5, u«y] - J>[w, u]) 

= — — / (x — V<^(x)) • C° V^(x) «o(x) dx + — / w(x) (w(x) — t'o(x)) dx 
hX JR d h J R d 

/ Tr(VC(x)) u m (x) dx - / u(x)((x) ■ Vv(x) dx 

X JR d JR d 

u(x)w(x)dx+ / [— Av(x) + av(x)] w(x) dx . 

JR d 

The above inequality being valid for arbitrary (C,w) & C^ D (M d ;M d ) x C£°(R d ), it is 
also valid for (— C, —w) so that we end up with 

(56) - / C ■ ( V (u m ) - x u Vv) dx + / It v ~ v ° _ At> + av - u ) w dx 

X JR d JR d \ h J 

1 f 

= 1 — / (x — V</?(x)) ■ C°Vv?(x) u (x) dx . 
hXJR d 

Observe that, since V (u m ) = 2 u m l 2 V {u m > 2 ) and u e L m (M d ) with V (u m / 2 ) e 
L 2 (K. d ) by the definition of /C and Proposition [8], Vu m belongs to L 1 (R d ) and the first 
term on the left-hand side of fl56l) is meaningful. 

Now, taking £ = in fl56l) and using a density argument, the regularity of the 
minimisers, see Proposition [8] and Corollary [9J readily gives ( l50l) . We next take 



IS 
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w = in fl56l) to obtain 



(57) / (■ (V(w m ) -x u Vv) dx = - / (x - Vtp(x)) ■ C° Vy»(x) ito(s) dx 

for allCeC oo (M a! ;M a! ). 

On the one hand, the Cauchy-Schwarz inequality and the properties of Vy? ensure 
that 

\ 1/2 

|C(x)| 2 u(x) dx J 



(x — Vyj(j)) • C° V<^(x) Mo(^) dx < W^w, u ) 
and we infer from the Holder and Sobolev inequalities that 
(x — Vcp(x)) ■ C° Vy?(x) mq(x) dx 



< W 2 (u,U )||C|| 2m Ih" 172 !!^-!) 

< CW 2 (« )M ,) nci| 2m ||v(W 2 )||^ /m . 

Recalling (157)) we arrive at 

/ C-(V(« m )-x«Vi;) dx <CW 2 ( Ml «,) HCIU \\V {u m/2 )\\ l 2 /m 

for all C G Cg°(R d ;]R d ), whence V (w m ) - x«Vt> G L 2 ™^ 2 ™- 1 ) (R<*) and (jSTJ) by a 
duality argument. 

On the other hand, consider £ G C^°(lR d ). Using Taylor's expansion, we have 

|£(x) - t(y<p{x)) - (V£ o V^)(x) • (x - cp(x))\ < Halloo |3; ~ V / (X)|2 

for x G M d . Multiplying the above inequality by u , integrating over M. d , and using 
the properties of V<p give 



[£n -£n-(V£oV^)-(id-^)n ] dx < p 2 £|U W 2 2 (u, u ) . 
Combining the above inequality with f )57|) (with ( = V£) leads us to fT^9|) . 

3. Convergence 



D 



Let x G (0, x c ), (uoj^o) G /C, and /i G (0,1). We define a sequence (uh, n , Vh,n)n>o 
in /C as follows: 



(58) 

and, for each n > 0, 



(wft,,o,^,o) := (u ,v ) 



{uh,n+i, v h,n+i) is a minimiser in /C of the functional 



F h>n [u,v\ := — 



h T \\V ~ V h>n \\ 2 

X 



+ £ a [u,v] , (u, v) G K,. 



Recall that the existence of («/i, n +i,^/i,n+i) * s guaranteed by Proposition [7] since x G 
(0, Xc)- We next define a piecewise constant time dependent pair of functions (uh, Vh) 
by 

{u h ,v h ){t) := (u hjn ,v h>n ), te[nh,(n + l)h), n>0. 
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3.1. Compactness. By the analysis performed in Section [2], (uh(t),Vh(t)) belongs 
to K, for all t > and (uh, Vh) is endowed with several interesting properties which 
we collect in the next lemma. 

Lemma 11 (Uniform estimates). There is a positive constant C§ depending only on 
d, x, t~, a, u , and v such that, for < s < t, 



(i) 



£ a [u h (t),v h (t)] < £ a [u h (s),v h (s)] 



> t h T \\v h>n+ i - v h!n \\ 2 ) < 2£ a [u , v \ h 



,n=0 



X 



111 



U h (t)\\ m +\\Vv h (t)\\2<G 



[IV) 

and, for t > h, 



K(*)lll + / \x\ 2 u h (t,x)dx<C 5 (1+t), 



(v) 



h 



V (u™ /2 ) (s)||* + || (Av h - av h + u h )(s)j \\l ds < C 5 (1 + 1) 



vi 



(V(«T) - X« fc Vt; fc ) Wll^fct- 1 !) ds < Cb (i + *) ■ 



Proof. • Energy and moment estimates: Let n > 0. It follows from the properties 
of the minimisers of J 7 ^ in /C, see Proposition [3, that 



(59) 



2/? 



W|(Mfe, n+ l,Mfe, r 



T ||W/i, n +l - Vh,n\\2 



+4p(i,n+l,V+l] — 4Kni"M 



Consider s > 0, t E (s, oo), and set iV := [£//i] and v := [s//i]. If iV > 1 + u, we 
sum up ( 1591) from n = u to n = N — 1 to obtain 



2V-1 r 






Wf(«fc, f n-i,Ufc,„) ,. ||2 
h 7" ||Vft,n+l ~ v h,nh 



+ £a[u h ,N, V h)N ] < £ a [Uh,v, Vh,v] 



Since both terms of the left-hand side of the above inequality are non-negative by 
Lemma [5] and Lemma EJ we deduce (EJ) as 

£ a [u h (t),v h (t)) = £ a [uh,N,v htN ] < £a[u h ,u,v h ,„] = £ a [u h (s),v h (s)] , 

and (JTTJ) by taking s = [y = 0) and t — > oo in 



2V-1 

E 



^2 MM+ 1 ' U M i II II 2 
\~ T \\V hn+ i - V hn \\ 2 

X 



< 2h£ a [u h>v ,v h)V ] 
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The estimate ( 1ml) directly follows from the lower bound on S a , see Lemma [5] and 
Lemma EJ while we deduce (ITvl) from (In]) thanks to the inequalities 

JV-l 



>V 2 (5 , U h:N ) < yV 2 (S , U ) + ^ W 2 (Uh,n+l, W/i,n) 

n=0 

/JV-l 
< W 2 (5 , U ) + V^ [J2 W|(«h,n+1, «h,n 



,n=0 



and 



JV-l 



|^h,Jv||2 < |ko||2 + ^ Ikh.n+l - Wh,n||2 



n=0 



< 1 1 1?0 1 1 2 + ViV ( J^ ||Uh,n+l - «h,n II 2 
\n=0 

• Additional estimates: Assume now that t > h. By the properties of the min- 
imisers of J-"^ in /C, see Proposition El 

H V ( M M 2 +l)ll2 + h \\&V hjn+1 - av h . n+ i + Mfc,„+i||2 

2 

< - (H[u h ,n] - Ti[Uh,n+l]) +C 2 h (£ a [Uh,n, Vh,n] + £ a [Uh,n, Vh,n]' L/(m_1) ) 

+ r (||V^ !n ||2 + a \\v h>n \\l - ||V^, n +i||2 - a IK,n+i|| 2 ) ■ 
Summing from n — to N — 1, it follows from ( 134)) . (158|) . (jl]), and ( 1rv|) that 
, zv-i zv-i 

— Yl II V ( M M 2 +l)ll2 + h Yl W Av h,n+l - (XV h , n+1 + U h>n+1 \\l 
* n =0 n=0 

2 7V_1 

< - (H[u h ,o] - H[u h , N }) + C 2 hJ2 (£a[Uh,n, Vh,n) + £ a [u h , n , V h ^ m '^) 
X n=0 

+ r (||Vvfc i0 ||2 + oi \\vh,o\\l ~ ||Vvfc,jv||| - a ||v/i,jv||1) 
<- [H[u hfi ] + C 3 + \x\ 2 u htN (x) dx) +t (\\Vv ht0 \\l + a\\v hfi \\l) 

JV-l 

+ C 2 /l ^ (£ a [u M , Vh,o] + ^K,0, ^,o] 1/(m_1) ) 
n=0 

<C(l + t), 
whence (jvj). 
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Finally, setting -d := 2m/(m+ 1) G (1, 2), we infer from (15T1) . On]), (jvj), and Holder's 
inequality that 



7V-1 



^ Yl II V ( n M+i) "PM+i V V+i 



n=0 



\° 
\2m/(2m-l) 



JV-1 



< C| /i 1 ^ ^ W. 



2 l M /i,n+l) u h,n) 



n=0 



'JV-1 



v U 



m/2 
'h,n+l 

/N-l 



■Q/m 



H^l 



<C h 1 ^ I ^W%(u h>n+1 ,U h , nJ 



<c 



n=0 
(N+l)h 



£ I v {< 



m/2 
+1 



,n=0 
(2-tf)/2 



(2-*)/2 



v(% /2 (^: 



ds 



□ 



<C(l + t), 

which completes the proof. 

The following corollary guarantees the compactness of (uh), ( v h), and (Vt>/i) with 
respect to the space variable x: 

Corollary 12. There is a positive constant Cq depending only on d, \, r , a , u o, an d 
vq such that, for t > h, 

f (\\Ms)\\l m + \\^u h (s)\\l) ds < C7 6 (l+t), 



/ |M3)||^.d5 < C 6 (l+t) 
Jh 



Proof. We infer from the regularity of the minimisers of J-h,n i n ^i see Proposition [HI 
and the estimates of Lemma [TT] that 

-t rt 9 

m/2x 



(60) / \\u h (s)\\r ds < C / \\u h (s)\C ^(u71(s) ds<C(l+t). 

Jh Jh 2 

Next, since m G (1, 2), we have Vw^ = 2u h ~ V(m™ )/m and we infer from the 
estimates (jml) and (jvj) of Lemma [TT] that 



\Vu h (s)\\ 2 m ds< — 



m^ 



m/2N 



v(%1( s 



K(s)||^ m 



ds<C (1 + i). 



Finally, since m > 1, we combine the estimates (TTITT) — (rvj) of Lemma [TT1 with ([60]) to 
conclude. □ 

We now turn to the compactness with respect to time. 

Lemma 13 (Time equicontinuity) . There is a positive constant C-j depending only 
on d, x, t, a, u , and v such that, for < s < t, 

\\v h (t) - v h (s)\\ 2 < C 7 (yt^ + Vh), 

Wit) - Uh (s)\\ H - (d+2) < c 7 (i + tY m+1 ^ 2m (t- s + hf m -^ 2m . 



■n 
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Proof. Consider s > 0, t > s and set N = [t/h] and v = [s/h]. Either N — v and 
(uh,Vh)(t) = (uh,Vh)(s) which readily implies the result. Or iV > v and, on the one 
hand, we infer from the estimate dUl) of Lemma Qj] that 



N-l 



\v h {t) - v h (s)\\ 2 = \\v htN - Vh,uh ^ ^2 H u M+i ~ v h,nh 

n=u 

1/2 



'AT-l 



< VN - V ^ IK,n+l - v h,n\ 



J 2(N-u)hS a [u ,v ] ^ c iVT ^ +VJl 



T 



On the other hand, it follows from the Euler-Lagrange equation (|49p that, for £ e 
nd v <n < N - 1, 

£,{u h ,n+i ~ Uh,n) dx 



C^(R d ) and v <n < N- 1 



< /i / |Vf | |V (u™ n+1 ) - X Mfe,n+i Vf M+ i| dx 

W2(«fc,n+l>«h,n) 

+ IKIIw 2 .- 2 ' 

whence, after summing up these inequalities from n = u to n = N— 1: 



£(«/,(£) -«h(s)) dx 



AT-l 



C( w ?i,iv - w h)2/ ) dx 

£(Uh,n+l - Uh,n) dx 



n=v 
N-l 



< h ^2 H V ^ll 2m II V ( M M+l) - X U h , n+ i Vv ft ,n+l| 



2m/(2m-l) 



AT-l 



+ ||£||wa.« 5Z W 2 (Uh,n+h Uh,n) ■ 

Using the estimates (Imj) and (TvTj) of Lemma [IT], we obtain 

AT-l 

^ W|(ttfe,n+1, ttfe,n) <Ch, 



and 



JV-l 



^ 5Z II V ( U M+l) ~~ X M M+1 V^,n+l| 



2m/(2m-l) 



AN+l)h 

< I II V «) -X% V^|| 2W(2rrl _ 1) ds 

J(u+l)h 

< C [1 + (TV + i)/ l ](-+ 1 )/ 2m [(TV - z/)^]^- 1 )/ 2771 

< C (1 + t)( m+1 )/ 2m (i - s + h) (m - 1)/2m , 
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so that 

_ £K(t) - u h (s)) dx < C (1 + tf m+ ^ 2m (t-s + hf m -^ 2m \\V£\\ 2m 

+ C ll^llw 2 .- h. 

Since H^R 01 ) and H d / 2 (IR d ) are continuously embedded in L°°(IR d ) and L 2m (M d ), re- 
spectively, we end up with 



£(u h (t) -u h (s)) dx 



< c (i + t f m+i)/2m (t- s + h) {m - i)/2m yen 



H d+2 



from which the result follows by a duality argument. □ 

3.2. Convergence. We are now in a position to establish the strong convergence of 
(uh)h and {vh)h required to identify the equations solved by their cluster points as 
h — > 0. For further use, we define the weight g by 

g(x) := — = , xeR d . 

Proposition 14 (Convergence). There are a sequence (hj)j>i of positive numbers 
in (0, 1), hj — > 0, and functions 

u E C([0, oo); LT (d+2) (R d )) , v E C([0, oo); L 2 (M d ; e (x) dx)) 

such that, for allt>0,5<E (0, 1), and T > 1, 

u hj {t) — >u(t) in rH d+2 )(M d ), 

Vhj{t) — ► v(t) m L 2 (M d ; p(x) dx) , 

u h .^u in L 2 (5,T;L m (M d )), 
(61) 

v hj — y v in L 2 (5,T; H 1 ^; g(x) dx)) . 

Proof. Since H 1 (M a! ) is continuously embedded in L 2 (M d ; g(x) dx), it follows from the 
estimates (lull) and flrvl) of Lemma [TTl that {t'h(^) : {t,h) £ [0, oo) x (0,1)} lies in a 
compact subset of L 2 (M d ; g(x) dx) while the time equicontinuity of Lemma [TBI entails 
that 

limsup / M«,x)- w( ,,x)|V») d,<C 7 VF^, 0< S <*. 

We then infer from a refined version of the Ascoli-Arzela theorem, see [H Proposi- 
tion 3.3.1], that there are a sequence (hj)j>i of positive numbers in (0, 1), hj — > 0, 
and v e C([0, oo); L 2 (M d ; p(x) dx)) such that, for all t > 0, 

v h .(t)-^v(t) in L 2 (R d ; g(x) dx) . 

It next readily follows from the estimate (Iiv|) of Lemma [Til and Lebesgue's dominated 
convergence theorem that 

(62) v hj — >v in L p (0,T;L 2 (R d ;g(x)dx)) 

for all p 6 [1, oo) and T > 0. 
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Furthermore, recall that H 2 (M d ) is compactly embedded in H 1 (R d ; g(x) dx) which 
is continuously embedded in L 2 (M a! ; g(x) dx). Given 5 G (0,1) and T > 1, {vh)h 
is bounded in L 2 (<5, T; H 2 (M d )) by Corollary [T2l We may then apply a standard 
compactness result, see [2H Lemma 9], to deduce from (1621) the second convergence 
stated in flBTT) . 

Next, H 1 / m (M d ) is densely and continuously embedded in L m ^ m ~ 1 \M d ). Conse- 
quently, L m (M d ) is continuously embedded in H -1 / m (IR d ) and we then conclude that 
L m (M d ) n L 1 ^; (1 + \x\ 2 ) dx) is compactly embedded in B.~ {d+2) iW 1 ) . Owing to the 
estimates ( fmj) and (trvT) of Lemma HU there is for each T > a compact subset of 
H" (d+2) (M d ) (depending on T) in which u h {t) lies for all h G (0,1) and t G [0,T]. 
This fact along with the time equicontinuity of (uh) in Lemma d3] and [U Propo- 
sition 3.3.1], ensure that there are a subsequence of {hj)j>\ (not relabelled) and a 
function u G C([0, oo); fH d+2 )(M d )) such that, for all t > 0, 

(63) u hj (t)-^u(t) in H"( d+2 )(M d ). 

Thanks to the estimate (Jml) of Lemma [TT] and the continuous embedding of L m (IR d ) 
in H^ d+2 )(IR d ), we may combine (1631) and Lebesgue's dominated convergence theorem 
to deduce that 

(64) u hj —>u in L p (0,T;H- (d+2) (M d )) 

for all p G [1, oo) and T > 0. 

Finally, given 5 G (0, 1) and T > 0, it follows from the estimates (Jml) -(Jiv1) of 
Lemma [11] and Corollary [12] that 

(65) (u hj )j>i is bounded in L 2 (5,T; W 1,m (R d )) n L 1 ^; (1 + |a;| 2 ) dx)) . 

Since W 1 ' m (M d )nL 1 (M d ; (1 + \x\ 2 ) dx) is compactly embedded in L m (R d ) and L m (M d ) 
is continuously embedded in H~^ +2 ^(]R d ), a further use of [2U Lemma 9] allows us 
to deduce from (|64"1) and (J65l) the first convergence stated in (|6"T1) . D 

3.3. Identifying the limit. In addition to the convergence of Proposition [14], we 
note that the estimates of Corollary [T2l ensure that, after possibly extracting a further 
subsequence, we may assume that 

(66) Av hj -^Av in L 2 ((<5,T) x R d ) , 

(67) u/y-^u in L 2m {5,T;L 2 (R d )). 

Now, checking that (it, v) is a weak solution to (CD) follows from Lemma dO] Propo- 
sition [TH (j66l) . and (J671) by classical arguments. Observe in particular that f lBTT) 
and ( ISTj) imply that (a^V^.) >1 converges weakly towards wVf in L 1 ((5, T) x 

R d ; £>(x) dx) for all 5 G (0, 1) and T > 1. We have thus constructed a weak so- 
lution to the parabolic-parabolic Keller-Segel system (J2J) and proved Theorem [TJ 
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